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Abstract. We present XMM-Newton observations of the O supergiant £ Orionis (09.7 lb). The spectra are rich 
' in emission lines over a wide range of ionization stages. The RGS spectra show for the first time lines of low ion 

stages such as Cvi, Nvi, Nvn, and Ovn. The line profiles are symmetric and rather broad (FWHM w 1500 km 
s" 1 ) and show only a slight blue shift. With the XMM-EPIC spectrometer several high ions are detected in this 
star for the first time including Arxvn and Sxv. 
Simultaneous multi-temperature fits and DEM-modeling were applied to the RGS and epic spectra to obtain 
' emission measures, elemental abundances and plasma temperatures. The calculations show temperatures in the 

, range ~0.07 - 0.6 keV. According to the derived models the intrinsic source X-ray luminosity at a distance of 251 

pc is i x = 1.37(.03) x 10 32 ergs" 1 in the energy range 0.3 - 10 keV. In the best multi-temperature model fit, 
the abundances of C, N, O, and Fe are near their solar values, while the abundances of Ne, Mg, and Si appear 
somewhat enhanced. 

The sensitivity of the He-like forbidden and intercombination lines to £ Ori's strong radiation field is used to 
qq ' derive the radial distances at which lines are formed. Values of 34 i?» for Nvi, 12.5 R* for Ovn, 4.8 R, for Neix, 

f^*) . and 3.9 R* for Mgxi are obtained. 
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1. Introduction hetg data Waldron & Cassinelli (2001) concluded that 

most of the emission is at energies below 1 keV. 

The object C Orionis is a supergiant of spectral type ft ig generally believed that most of these x _ rays 

09.7 lb. Its X-ray spectrum has previously been ob- TecfAved from stars origiliat e in shocks distributed 

served with various satellites such as Emstem (Cassinelli throughout thc wind; and that these shocks are created 

& Swank 1983), ROSAT (Berghofer et al. 1997), and by inst abilitie S in the line-driving mechanism which drives 

ASCA (Kitamoto et al. 2000). The Emstem observations ^ QveraU outflow , Lucy & Solomon 1970; Lucy 1982; 

were particularly interesting: the Solid-State Spectrometer 0wocki et ^ ^ Degsart & Qwocki 20Q5) Tq explore 

observations of C Ori were the first to show hot star X- &e relatiQn between ^ Xray emission from ( Qri &nd 

ray emission lines, those of Sixm and Sxv (Cassinelli & ^ ties of the shocks in its wind line files from 

w 11 ' ' high-resolution Chandra observations were used by Cohen 

From R OSA T observations Berghofer et al. (1996) e t al. (2006), Owocki & Cohen (2006), and Oskinova et al. 

found that the X-ray emitting plasma could be character- (2006) 

ized with an overall kT of 0.22 keV and that the X-ray to Waldron & Cassinelli (2001) also analyzed the ra- 
Bolometric luminosty ratio of log L K /L hol = -6.74 is typi- tios of thc f orbidden (/) 5 the intercombination(z) and thc 
cal for O stars and early B supergiants. Based on Chandra- resonance ( r ) emission lines in the He-like ions of Si XIII, 
Mgxi, Neix, and Ovn. The lines were found to be pro- 
Send offprint requests to: A.J.J. Raassen e-mail: duced at a wide range of radii. The specific radius for each 
a. j . j .raassen@sron.nl He- like ion was found to correspond to X-ray optical depth 
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Table 1. Observation log of the data of C Ori. 



Instr. 


Filter 


Mode 


Date-obs-start 


Date-obs-end 


Duration(s) 


MOS1 


Thick 


Small Window 


2002-09-15T13:12:32 


2002-09-16T00:51:00 


41728 


MOS2 


Thick 


Timing 


2002-09-15T13:12:28 


2002-09-16T00:46:41 


41473 


pn 


Thick 


Full Frame 


2002-09-15T14:04:08 


2002-09-16T00:51:20 


38362 


RGS1 


None 


Spec+Q 


2002-09-15T13:ll:17 


2002-09-16T00:52:38 


41979 


RGS2 


None 


Spec+Q 


2002-09-15T13:ll:17 


2002-09-16T00:52:38 


41979 



o 




1 O 



20 

Wavelength (A) 



30 



Fig. 1. Background-subtracted epic-mos and first-order RGS spectra of C Ori. Error bars indicate la statistical errors 
including the effect of background. The grey (red in the electronic version) curve shows the best-fit model. A number 
of prominent lines are labeled with the emitting ions. 



unity to the continuum opacity at each line wavelength. 
Their most surprising result was that the high ion stages, 
such as Si xni, appear to form very close to the star. At 
these close distances the shock velocity jump was found 
to be larger than could be explained by the speed that the 
wind should have at these small radii. 

Recently Pollock (2007a) developed a new paradigm 
for O-star X-ray emission, which states that charge ex- 
change, ionization, and excitation are likely to be pro- 
duced by protons far out in the wind. However, this ap- 
proach is not yet capable of describing the spectrum and 
line fluxes in a quantitative way, and so it will not be 
discussed further in this paper. 

Variability has been found in both ( Ori's Ha profile 
(Ebbets 1982; Kaper et al. 1997) and its UV P-Cygni pro- 
files (Snow 1977; Kaper et al. 1996; Kaper et al. 1999). 
However, the X-ray light curve has provided no evidence 
for variability over a long time range (Berghofer & Schmitt 
1994). 

In many of the recent papers the authors have focussed 
on line profiles and individual line features. This is espe- 
cially true for spectra obtained by means of the high res- 
olution instruments aboard Chandra (hetg and meg). In 



contrast RGS aboard XMM-Newton provides us with new 
information regarding the longer wavelength part of the 
spectrum where ion stages such as Cvi, Nvi, and Nvn 
can be observed. In addition, the higher collecting area 
of XMM-Newton allows us to obtain a better estimate 
of the distribution of emission measure versus tempera- 
ture. In the present paper the spectra and line-fluxes of 
£ Ori, observed with XMM-epic and -RGS, are described 
for the first time by means of a multi-temperature fit and 
DEM-modeling to the total spectrum. The approach we 
choose to take is to use emissivity results from optically 
thin Coronal Ionization Equilibrium (CIE) plasma calcu- 
lations. 

Apart from this broad description of the total spec- 
trum, individual line fluxes and line ratios are analyzed in 
this paper. Ratios of the fir lines of He-like ions are used 
to determine the typical distances between the X-ray emis- 
sion regions and the stellar photosphere. The wavelength 
sensitivity of XMM-Newton allows us to check the results 
of Waldron & Cassinelli (2001) for Mgxi, Neix, and Ovn 
(with somewhat more precise radial ranges) , while extend- 
ing the analysis to include the Nvi lines. The results for 
the lines of N vi are particularly interesting because this 
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is the first time that this line ratio has been interpreted 
in this manner for this star. 

( Ori was observed on 2002 September 9 by means 
of XMM-Newton. A log of these observations is shown 
in Table 1. The spectral analysis by means of multi- 
temperature fitting and DEM-modeling are described in 
Sect. 2. The analysis of the line profiles and the individual 
line fluxes, including several that have not been studied 
earlier, is given in Sect. 3. The results are discussed in 
Sect. 4. Parameters of £ Ori are given in Table A.l in the 
appendix. 

2. Spectral analysis 

2.1. Multi-temperature fitting 

We have determined the thermal structure and the el- 
emental composition of £ Ori's X-ray emitting plasma 
by simultaneously fitting multi-temperature models to 
the RGS, epic-mosI, and pn spectra. The fit was per- 
formed using our spectral analysis program SPEX (Kaastra 
et al. 1996a) in combination with the MEKAL (Mewe- 
Kaastra-Liedahl) code, as developed by Mewe et al. (1985, 
1995) (see Fig. []}. mekal calculates a continuum and 
models more than 5400 spectral lines. It is available on 
the wwv\0. The applied model is a Collisional Ionization 
Equilibrium (CIE) model for optically thin plasma. The 
ionization equilibrium is based on calculations by Arnaud 
& Rothenflug (1985) and Arnaud & Raymond (1992). 

2.1.1. Temperature and Emission measure 

In the multi-temperature calculations we used initially five 
temperatures. Out of the five fitted temperature compo- 
nents three turn out to be significant. These three tem- 
perature components were coupled to one Nn absorption 
column density and one set of abundances, which were 
free to vary. The column density, temperatures, emission 
measures EM = J n e nndV , X-ray luminosities, elemental 
abundances, and micro-turbulent velocity (v m ic) are given 
in Table 2, together with the statistical ltr uncertain- 
ties. Three temperatures at T=0.073 keV, T=0.201 keV 
and T=0.551 keV were determined. Interpreting the re- 
sults in terms of a wind with embedded shocks and using 
kT = 1.2 [u(1000 kms -1 )] 2 keV, we find that the kT 
range from 0.073 keV to 0.551 keV corresponds to shock 
jumps from 246 km s _1 to 676 km s~ x . These are well be- 
low the terminal wind speed of 1885 km s _1 (= 2.3 keV), 
so they seem physically reasonable. 

Out of the three EM values the most uncertain is the 
one corresponding to the lowest temperature component 
(T=0.07 keV). This component is determined by the C vi 
and Nvi features. The most robust of the three compo- 
nents is the one at 0.20 keV because it involves the largest 
number of emission lines in the spectrum. The EM values 
for £ Ori decrease with T and approximately follow the 
relation EM cx T'°- s . 



Table 2. Multi-temperature fitting for combined epic-pn, 
epic-mosI & RGS spectra of £ Ori a . 



comp. i 


value 


Nh" 


5.0 




feTi c 


0.551±0.013 




EM\ d 


1 57±0 18 






4.18 






3.8 




Umicl 9 


1080±135 




kT 2 c 


0.201±0.004 




EM2 d 


3.76±0.31 




ix2 E 


6.71 






20.6 




«mic2 9 


950±70 




kT 3 c 


0.073±0.006 




EM 3 d 


10.1±2.7 






2.8 




w 


111.0 






1240±230 




J2i EM i d 


15.4±2.2 






13.7±0.3 






135±4 






-3 10" 4 






Abundance 8 






Ai 


Ai/Ao 


C 


1.02±0.11 


1.20±0.05 


N 


1.03±0.09 


1.21±0.05 


O 


0.85±0.07 


1 


Ne 


1.33±0.15 


1.57±0.07 


Mg 


1.66±0.16 


1.96±0.08 


Si 


1.30±0.14 


1.54±0.05 


Fe 


0.99±0.09 


1.15±0.06 


Xred 


1246/1064=1.17 



Notes: 

a: HIPPARCOS distance d = 0.25 kpc is used. 

b: Column density (10 20 cm -2 ) derived for standard 

absorption model (with solar abundances) from Morrison 

& McCammon (1983) for all components, 
c: Temperature in keV. 
d: Emission measure in 10 54 cm -3 . 

e: X-ray luminosity of the instrumental XMM-Newton band 

(0.3-10 keV) in 10 31 ergs -1 . 
/: X-ray luminosity of the total emission range, 

0.03-0.9 keV (see Fig. HI, in 10 31 ergs" 1 . 
g: Micro-turbulent velocity parameter used to characterize 

line broadening (in kms" 1 ). 
h: Artificial redshift to match the observed wavelengths to 

the theoretical ones. 
i: Abundance relative to solar photospheric number 

abundance (Anders & Grevesse 1989 for all ions except 

Fe, where Grevesse & Sauval 1998, 1999 was used). 



1 http:/ /www. sron.nl/divisions/hea/spex/ 
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The total EM over the three temperature bins is 1.35 x 
10 55 cm~ 3 . It is interesting to compare these values to 
the "wind emission measure" EM W given by Cassinelli et 
al. (1981). Based on their formula we obtain for £ Ori 
EM W = 2.44 x 10 58 cm~ 3 . This shows that only a small 
fraction of the wind is contributing to the X-ray spectrum. 

The emission lines observed in £ Ori's X-ray spec- 
trum clearly exhibit Doppler broadening of roughly 1000 
± 100 kms -1 . Therefore, the lines generated in the model 
spectrum need to be broadened in excess of instrumen- 
tal effects to allow a reasonable comparison with the 
data. For this purpose, we used the micro-turbulent ve- 
locity parameter (v m i C ) available in SPEX. In using this 
parameter, we do not mean to imply that the emission 
line broadening is microturbulent in nature, rather, it is 
just a convenient way to parameterize the line broad- 
ening due to the outflow of the X-ray emitting regions. 
The v m i C parameter is related to the FWHM measure of 
the line broadening through the relation FWHM/E t, a = 
2\/ln 2v m i c /c. Here E D b s is the energy of the emission line, 
i.e. 12.39842/Aobs. A characteristic value used in the fit- 
ting was FWHM/E obs = 5.5 ± 0.5 x 10" 3 . 

To achieve the best fit possible, the overall Doppler 
shift parameter z rec i was used. The best fits were achieved 
using an overall blueshift of roughly 100 kms -1 . To inter- 
pret this number, it is interesting to note that some small 
blueshift of the spectral line centroids is to be expected 
due to the attenuation of the X-rays emitted from shocks 
moving away from the observer on the far side of the wind. 
Nonetheless, the fact that this shift is roughly of the order 
of systematic uncertainties in the wavelength calibration 
(Pollock 2007b) means this result should be viewed with 
some caution. The interpretation of individual emission 
line shapes is considered in more detail in Sect. 3. 

2.1.2. Abundances 

The elemental abundances are relative to solar photo- 
spheric values from optical studies (Anders & Grevesse, 
1989) except for Fe, for which we use log Ap c = 7. 5C0 (see 
Grevesse & Sauval 1998 and 1999) instead of 7.67 (Anders 
& Grevesse 1989). 

The abundances have been determined through a fit to 
the total spectrum, allowing them to freely vary along with 
the emission measures and temperatures. When interpret- 
ing these abundance measurements, two issues should be 
kept in consideration: First, the abundances are strongly 
anti-correlated with the emission measures (especially 
in cases involving a weak continuum), with the result 
that the products of emission measure and abundance 
(EM x Ai) are quite robust. Second, the ratios between 
the abundances (Ai/Ao or Ai/Ap c ) of the various elements 
are more robust than the absolute abundances themselves. 
Therefore, the ratios between the abundances are often 
presented in papers. Here we give both values, the abso- 

2 Here log Ap e is the logarithm of the Fe-abundance relative 
to log^ H = 12.0. 



lute abundances and the ratio of the elemental abundance 
over the oxygen abundance averaged for a number of fits. 

The carbon abundance is strongly related to the low- 
est temperature and the corresponding emission measure. 
Almost 50% of the line flux of this ion is produced by 
the lowest temperature component. The obtained abun- 
dances are all close to solar photospheric values except for 
Ne, Mg, and Si. This behaviour is different from coronal 
plasmas, for which in the quiescent state a First Ionization 
Potential (FIP) effect is noticed. This effect implies that 
elements with low first ionization potential (Mg, Si, Fe) 
are overabundant in the corona compared to the photo- 
sphere. For active, flaring coronal plasma an inverse FIP 
effect is often found. 

2.1.3. X-ray luminosity 

The luminosity is a robust quantity. The X-ray luminosi- 
ties, given in Table 2, are the model luminosity in the 
band 0.3 - 10 keV (the energy range of XMM- Newton) 
and in the energy band 0.03 - 0.9 keV, the energy range of 
the determined emission measure (see Fig. [2]). Both val- 
ues correspond to the location of the emitting plasma, 
i.e., they are corrected for absorption by the interstellar 
medium (ISM) and by the stellar wind. Our instrumental 
(0.3 - 10 keV) value of L x = 1.37 x 10 32 ergs" 1 results in 
a value for logL x / L\, i = -6.66(.09) which is similar to the 
values obtained using Einstein (Chlebowski et al. 1989) 
and ROSAT (Berghofer et al. 1996). 

In order to attribute the X-ray emission from this star 
to shocks embedded in its wind, some of the mechanical 
energy of the wind outflow must be converted into heat 
energy. It is therefore interesting to compare the mea- 
sured X-ray luminosity to the mechanical luminosity of 
the wind, a quantity which can easily be calculated from a 
formula given in Howk et al. (2000). The resulting value of 
L w — 1.6 x 10 36 erg s _1 should be compared with the X-ray 
luminosity of the total emission of L x = 1.4 x 10 33 erg s _1 . 
This indicates that only about a one thousandth of the 
mechanical energy available in the wind is radiated away 
as X-rays, a small fraction. 

2.2. DEM modeling 

Apart from a multi-T fitting, the combined EPIC-MOS, 
epic-pn, and RGS spectrum was also fitted by means 
of a differential emission measure (DEM) modeling (e.g., 
Kaastra et al. 1996b) using the regularization module in 
SPEX. This module constructs an emission measure distri- 
bution with the constraint that the first and second deriva- 
tives are smooth and continuous. We define the DEM by 
n e nndV/d\ogT (integrated over one logarithmic tempera- 
ture bin, the emission measure = n e n}iV). 

The result for the emission measure distribution EM 
is shown in Fig.[2j The temperature covers the range from 
0.03 keV to 0.9 keV, i.e., from 0.3 to 10 MK, with three 
temperature peaks around 0.06 keV, 0.2 keV, and 0.6 keV. 
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Fig. 2. DEM modeling of the epic-mos and RGS spec- 
tra of £ Ori. The emission measure EM per logarithmic 
temperature bin is in units of 10 58 cm~ 3 . 

These values are consistent with the results of the 3-T 
fittings. Just as in the case of the 3-T fitting the DEM- 
modcling is most robust for the range around 0.2 kcV, 
while the error bars are larger for the low temperature 
regime. 

3. Emission line fluxes 

As mentioned in the introduction, the shapes of the X- 
ray emission lines detected by high resolution X-ray spec- 
trometers have been very important in determining the 
nature of the X-ray emission from these stars. Therefore, 
apart from the description of the total X-ray spectrum of 
£ Ori by means of a CIE model, the fluxes, wavelengths, 
and broadenings of the individual emission lines were mea- 
sured. This was done by folding a Gaussian through the 
response matrix. 

In cases where line blending would interfere with a 
measurement, we adopted the following procedure to iso- 
late the contribution of the line being measured: The pa- 
rameters derived in the multi-T fit (Sect. 2.1.1.) were used 
to generate a model spectrum which included the contri- 
butions of all relevant ions except the ionization state of 
the line being measured. A single Gaussian (whose param- 
eters were free to vary in the fit) was then added on top of 
the computed model spectrum to measure the individual 
contribution of the line being measured. 

In Table 3 the wavelengths, the measured fluxes at 
Earth (not corrected for interstellar absorption), and the 
widths of a number of prominent lines are collected. 
Within brackets the la errors are given. The rest wave- 
lengths (Ao) are taken from Kelly (1987) and Dere et al. 
(2001), 

3.1. Analysis of emission line profiles 

The lines in the RGS spectrum of ( Ori are symmet- 
ric, broadened and slightly blue-shifted. After correcting 
for instrumental broadening, FWHM is s=a 1500 km s _1 



(cf. Tabic 3). For the dominant, unblended lines of Cvi 
(33.736 A), Nvn (24.781 A), Ovm (18.969 A), Fexvn 
(15.013 A), and Ncx (12.134 A) wc obtain a FWHM of 
1540(240) kins" 1 , 1500(360) kms" 1 , 1600(140) kms" 1 , 
1440(260) kins" 1 , and 1400(600) kms" 1 , respectively. 
Our results confirm the observations done with the 
Chandra-HETGS (Waldron & Cassinelli 2001, Miller 2002, 
Miller et al. 2002). The line broadenings of the individual 
line measurements are consistent with the relative broad- 
ening FWHM/E obs of 5.5 (in units 0.001), determined 
by means of the v m ; c value of the multi-temperature fit in 
Sect. 2.1.1. For the same lines as mentioned above the 
values of FWHM/E obs are 5.2(0.8), 5.0(1.2), 5.4(0.5), 
4.8(0.8), and 4.7(1.9), respectively. This broadening is con- 
sistent with a picture of outwardly moving shocks. 

Most wavelengths of the lines in Table 3 show a nega- 
tive deviation, corresponding to w -90(60) kms -1 , which 
is based on the same five lines as used above with indi- 
vidual deviations of -200(90) kms" 1 , -120(110) kms" 1 , 
-75(50) kms-\ -20(80) kms" 1 , and -70(200) kms" 1 . 
There is some tendency that the features of the lower ion- 
ized (cooler) ions are more blue-shifted. However, due to 
the large errors, the evidence for blue-shifted line profiles 
is not very strong. 

The possible effects of resonance line scattering opac- 
ity on the shapes of hot star X-ray emission lines have 
been discussed in Leutenegger et al. (2007) and Ignace & 
Gayley (2002). We use the measured line fluxes of FeXVll 
to estimate the optical thickness of the shock sources of 
£ Ori. In an optically thin plasma, the line flux ratio of the 
Fexvn lines at 15.013 A and 15.260 A would be expected 
to be 3.5(1.0) from the MEKAL atomic data, while the 
measured ratio for our ( Ori data is 2.5(0.5). Thus we can 
conclude that there is no severe resonance line scattering 
in the FeXVll line at 15.013 A and that the individual line 
source regions are optically thin, even though the entire 
wind is thick to X-rays at these wavelengths. 

The measured line fluxes of the He-like lines of Nvi, 
O VII, Neix, and Mgxi are used in the next section. 

It is likely the weak features in the EPIC-MOS spectrum 
near 3.9 A and 5.0 A are due to Arxvn and Sxv, respec- 
tively. The optimal formation temperatures of these ions 
(22 MK or 1.9 kcV, and 16 MK or 1.4 keV, respectively) 
are much higher than any of the temperatures found in 
our three-temperature and DEM temperature modeling 
in Sect. 2. The presence of emission lines from these ions 
indicates that they are formed in a high energy tail in the 
emission measure distribution of this star and in the wings 
of the line profiles. 

3.2. He-like line diagnostics 

Gabriel & Jordan (1969) were the first to demonstrate that 
forbidden (/), intercombination (i) and resonance (r) lines 
of the He-like "complexes" provide useful diagnostics for 
X-ray emitting plasmas. The most recent studies includ- 
ing the effects of dielectronic satellite lines and a radiation 



6 A.J.J.Raassen et al.: XMM observations of (" Ori 

Table 3. Measured line fluxes at Earth of ( Orionis. Statistical la errors in the last digits in parentheses. 



ion 


K 


^obs 


AA 


Av c 




flux 


FWHM d 


HWHM e 




r i i 

[A] 


r % 1 

[A] 


r % ^ 

[A] 


n —It 

[km s x ] 


[10 ld 
cm 


erg 

— 1 T 


[10^ 3 kcV] [km s" 1 ] 


EPIC 


















Ar xvn 


3.949 


3.92(.031) 






0.07( 


.03) 


- 


- 


Sxv 


5.039 


5.032(.018) 


~ 




0.18( 


.05) 


- 


- 


Si xin 


6.692 


6.684 






0.84( 


.42) 


- 


- 


Mg xi 


9.170 


9.177(.014) 






1.18( 


.09) 


- 


- 


RGS 


















Si xin 


6.692 


6.688(.066) 


— 




0.84( 


.60) 


- 


- 


Mg xn 


8.423 


8.405(.030) 






0.27( 


.13) 


- 


- 


Mg xi 


9.169 


9.169 7 


— 


— 


0.59( 


.15) 


7.5(7.5) 


800(800) 


Mg xi 


9.231 


9. 23^ 




— 


0.31( 


.14) 


7.5(7.5) 


800(800) 


Mg xi 


9.313 


9.313 / 




— 


0.33( 


.12) 


7.5(7.5) 


800(800) 


Nex 


10.240 


10.212(.020) 


-.028 


-820(590) 


0.29( 


.14) 


<11.0 


<1330 


Nex 


12.134 


12.131(.008) 


-.003 


-70(200) 


1.83( 


.26) 


4.8(1.9) 


700(300) 


Fe xvn 


12.264 


12.256(.023) 


-.008 


-200(600) 


0.61( 


.19) 


4.8(1.9) 


700(300) 


Ne ix 


13.447 


13.444(.006) 


-.003 


-70(130) 


3.35( 


.39) 


3.8(1.0) 


620(160) 


Ne ix 


13.553 


13.548(.009) 


-.005 


-110(200) 


2.26( 


.32) 


3.8(1.0) 


620(160) 


Ne ix 


13.700 


13.699(.017) 


-.001 


-20(400) 


0.66( 


.20) 


3.8(1.0) 


630(170) 


Fe xvni 


14.205 


14.204(.011) 


-.001 


-20(220) 


0.79( 


.18) 


3.7(3.0) 


640(515) 


Fe xvn 


15.013 


15.012(.004) 


-.001 


-20(80) 


5.33( 


.33) 


4.0(0.7) 


720(130) 


O vin 


15.176 


15.176 / 






0.99( 


.44) 


4.0(0.7) 


720(130) 


Fe xvn 


15.260 


15.261(.014) +.001 


+20(280) 


2.16( 


.42) 


4.0(0.7) 


720(130) 


Fe xvn 


15.449 


15.415(.026) 


-.034 


-700(500) 


0.54( 


.26) 


<8.3 


<1550 


O vin 


16.006 


16.017(.008) +.011 +210(150)' 1 


2.21( 


.29) 


4.9(1.2) 


950(230) 


Fe xvni 


16.078 












- 


- 


Fe xvn 


16.775 


16.771(.006) 


-.004 


-70(110) 


2.04( 


.22) 


2.0(1.1) 


410(220) 


Fe xvn 


17.051 


17.05P 






2.84( 


.52) 


3.7(0.9) 


760(190) 


Fe xvn 


17.100 


17.100^ 






2.34( 


.42) 


4.7(1.3) 


970(270) 


Ovn 


17.396 


17.386(.034) 


-.010 


-170(600) 


0.62( 


.23) 


- 


- 


Ovn 


17.768 


17.740(.014) 


-.028 


-500(250) 


0.64( 


.12) 


4.0(1.5) 


860(320) 


Ovn 


18.627 


18.617(.009) 


-.010 


-160(144) 


1.10( 


.16) 


3.9(0.8) 


880(180) 


O vin 


18.969 


18.962(.003) 


-.007 


-75(50) 


10.78( 


.36) 


3.5(0.3) 


800(70) 


Ovn 


21.602 


21.594(.005) 


-.008 


-110(70) 


8.07( 


.60) 


2.4(0.3) 


630(80) 


Ovn 


21.804 


21.775(.007) 


-.029 


-400(100) 


8.44( 


.63) 


3.8(0.4) 


1000(110) 


Ovn 


22.101 


22.088(.009) 


-.013 


-180(120) 


0.93( 


.15) 


3.0(0.7) 


800(190) 


N vn 


24.781 


24.771(.009) 


-.010 


-120(110) 


1.44( 


.14) 


2.5(0.6) 


750(180) 


Nvi 


24.898 


24.898(.041) 


-.000 


0(490) 


0.34( 


.16) 


1.8(0.6) 


530(180) 


Caxni 


26.719 


26.712(.055) 


-.007 


-80(640) 


0.13( 


.08) 


<2.4 


<790 


Cvi 


26.990 


26.972(.024) 


-.018 


-200(270) 


0.24( 


.09) 


<2.2 


<720 


Arxiv? 


27.41 


27.343(.029) 


-.067 


-700(320) 


0.19( 


.10) 


<5.0 


<1600 


Caxn 


27.973 


27.955(.027) 


-.018 


-190(300) 


0.25( 


.10) 


1.6(1.3) 


540(440) 


VI 


Zo.400 


28.396(.031) 


-.070 


-700(300) 




. 1Z 1 


2.9(1.5) 


1000(520) 


Nvi 


28.787 


28.771(.027) 


-.016 


-170(280) 


0.79( 


.12) 


3.0(0.6) 


1040(210) 


Nvi 


29.084 


29.076(.027) 


-.009 


-90(280) 


0.97( 


.14) 


3.0(0.6) 


1050(210) 


Nvi 


29.534 


29.521(.059) 


-.013 


-130(600) 


0.21( 


.07) 


3.0(0.6) 


1070(210) 


Sxiv 


30.441 


30.427(.023) 


-.014 


-140(230) 


0.16( 


.08) 


<1.0 


<370 


Si xn 


31.016 


30.959(.033) 


-.057 


-550(320) 


0.17( 


.07) 


<2.4 


<900 


Sxm 


32.22 


32.226(.030) 






0.22( 


.08) 


<1.7 


<660 


Sxiv 


32.554 


32.527(.043) 


-.027 


-250(400) 


0.28( 


.09) 


<2.2 


<870 


Cvi 


33.736 


33.713(.010) 


-.023 


-200(90) 


2.96( 


.33) 


1.9(0.3) 


770(120) 


Sxm 


35.66 


35.622(.021) 






0.53( 


.08) 


1.5(0.5) 


650(220) 


Sxn 


36.398 


36.350(.030) 


-.048 


-400(250) 


0.40( 


.10) 


1.3(0.7) 


570(310) 



Notes: 

a: Ao is theoretical wavelength from Kelly (1987) and Dere et al. (2001) and A b s is observed wavelength with the statistical 
la error in parentheses; b: line wavelength shift; c: velocity shift; d: line broadening fwhm; e: line half-width at half 
maximum; /: line wavelength fixed to assumed "blue-shifted" value; h: if we take the average wavelength for the Ovin, 
Fexvin blend then we obtain a velocity shift of -500 kms -1 ; i: line wavelength fixed to theoretical value. 
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Fig. 3. The panels show the He-like O vn and Ne ix line triplets consisting of a resonance line (r), an intercombination 
line (i), and a forbidden line (f) for the 09.51b star £ Ori. 

In both panels, the forbidden line nearly disappears due to radiative depopulation of the upper level in the strong 
radiation field relatively close to the surface of the star. 



field have been done by Porquet et al. (2001). The ratio 
f ji is dependent on electron density because of the de- 
population of the upper level of the / line in favour of the 
upper level of the i line at increasing density. This effect 
can be used to derive electron densities in circumstances 
where radiation fields are relatively weak, such as in cool 
stars. However, in OB stars, where the radiation fields are 
much stronger, the depopulation of the upper level of the 
/ line in favour of the upper level of the i line occurs by 
radiative absorption (e.g., Blumcnthal et al. 1972, Porquet 
et al. 2001). The f/i ratio no longer indicates the density, 
but instead provides information on the mean intensity 
of the radiation field, hence the radial distance R of the 
X-ray source from the star. This is illustrated in Fig.[3]for 
the Ne IX and O vn triplets. It is seen that the forbidden 
line is nearly suppressed for both ions. 

The fact that the UV radiation field is the dominant 
effect for £ Ori is confirmed by detailed calculations. Using 
the formalism developed by Blumenthal et al. (1972), we 
calculate the radial dependence of R in the envelope of 
£ Ori on the basis of a Planck curve (Ness et al. 2001) of 
31500 K as well as on the basis of a photosphcric model 
of Vacca et al. (1996) with T Gff =32000 K and log g=3.23 
and applying a UV flux model OSTAR2002 of Lanz & 
Hubeny (2003) with T Gff =32500 K, log g=3.25, and solar 
metallicity. These parameters closely correspond to those 
given in the appendix. The radial-dependent wind density 
is derived from the stellar and wind parameters of Lamers 
& Leitherer (1993) using a standard velocity law that is 
modified below 1.02i?„ (where i?* is the stellar radius) to 
provide a smooth density transition to the photospheric 
structure. The mean intensity of the UV radiation is large 
near the surface of the star and it decreases outwards by 

/ 9\ 1/2" 

l-(l-(^) 2 ) . Asare- 

suit of the radial dependence of the radiation field, the 
observed f/i ratio can be used to derive the radial loca- 
tion of the He-like ions that are producing the observed 
fir lines. For the four He-like ions which allow reliable 
measurements of their f/i line ratio (Nvi, Ovu, Neix, 



dilution factor W(r) = 



and Mgxi) these parameters indicate that ( Ori's radi- 
ation field will suppress the forbidden lines by radiative 
de-excitation of the upper level of / to much greater radii 
than would be possible with collisions (see, e.g., Waldron 
& Cassinelli 2001). Thus the expected f/i ratios as func- 
tions of radii for these ions are entirely controlled by the 
strength of £ Ori's UV radiation field. 

The predicted R dependencies of f/i are shown 
in Fig. [3J The f/i ratios derived from the RGS spec- 
tra for Nvi, Ovii, Neix, and Mgxi are: 0.22±0.12, 
0.114±0.027, 0.26±0.16, and 1.08±0.61, respectively (cf. 
Table 3) , agreeing with values for Ne ix and Mg xi given by 
Oskinova et al. (2006). These ratios correspond to average 
distances from the stellar surface R of the ion formation 
of 34±10i?*, 12.5il.5_R*, 4.8±1.8i?*, and 3.9±1.7i?*, re- 
spectively. For the last three ions this is in agreement with 
values established by Waldron & Cassinelli (2001), based 
on HETGS observations. They determined values for for- 
mation radii of <12i?„ for Ovu, of <^6R* for Neix, and 
3-5i?* for MgXI. 

The formation distances of He-like ions increase with 
decreasing ionization stage, but they increase also with 
increasing wavelength (A) of the line triplets. The forma- 
tion distance is proportional to some power of the wave- 
length of the intercombination line (A p , withp w 2.5—4.5). 
Waldron & Cassinelli (2001) pointed out that the derived 
radii (Rf ir ), at which the He-like lines are observed, are 
strongly correlated with the radii R\ at which the contin- 
uum adjacent to the line has an optical depth unity. Since 
the opacity varies roughly as A 3 this correlation would 
mean that we are observing the He-like lines from the 
deepest possible part of the wind from which X-rays could 
escape to the observer. 

Leutenegger et al. (2006) give estimates for minimum 
radii of formation of the hot plasma for O-stars. Their val- 
ues for minimum radii are about 1.5i?*, based on the ions 
of higher stages of ionization (Sxv, Sixni, and Mgxi). 
The reader should keep in mind that our measurements 
of the radii of line formation are not directly comparable 
to theirs. Our values are average values with a lower and 
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Fig. 4. The curves show the expected f/i line ratio as a 
function of formation radius R for the labeled He-like ions. 
The thickened portion of each line gives the radial range 
of the X-ray emission implied by the measured f/i line 
ratios. 

higher limit due to the uncertainty in the measured line 
fluxes and do not imply a minimum or maximum distance 
of ion formation. 

Finally, we used the resonance line, the forbidden line, 
and the intercombination line of the He-like lines as a tem- 
perature diagnostic for the individual regions responsible 
for the fir emission of each ion. As a single-temperature 
measurement method for specific regions in the wind, this 
method gives independent, complementary information to 
the global temperature fitting carried out in Sect. 2. Using 
the results of Porquet et al. (2001) we extract from the 
measured (i + f)/r ratios (cf. Table 3) for Nvi , Ovn , 
and Neix electron temperatures of about 0.05, 0.16, and 
0.24 keV, respectively. These three temperatures are in 
the range of the lower two of the three temperatures found 
with multi-temperature fitting in Sect. 2.1.1. 

4. Discussion of results and conclusion 

The O supergiant £ Orionis (09.7 lb) has a line rich X-ray 
spectrum. The spectral lines are symmetric and broad- 
ened. The FWHM corresponds to 1500 km s _1 ; less than 
the terminal speed of 1885 km s _1 derived by Lamers et 
al. (1999). The velocity widths are consistent with a pic- 
ture of outwardly propagating shocks that are distributed 
throughout the wind and that are produced by the insta- 
bilities associated with a line-driven wind. The best fits 
require a small blue-shift. It should be noted, however, 
that this blue-shift is of the order of the calibration un- 
certainties of the instrument. 

The continuum and line features are very well de- 
scribed by a 3-T collisional ionization equilibrium model 
(CIE) for optically thin plasma, using SPEX in combina- 
tion with MEKAL. The determined temperatures range 



from 0.07 to 0.55 keV (0.85 - 6.4 MK), corresponding to 
shock jumps from 246 km s _1 to 676 km s _1 . The total 
X-ray luminosity of the star corresponds to an emission 
measure that is far below the emission measure of the en- 
tire wind (Cassinelli et al. 1981, Cassinelli & Swank 1983), 
indicating that only a small fraction of the wind is partic- 
ipating in the X-ray emission. 

We used the values of the forbidden to intercombina- 
tion line ratio of He-like ions to derive characteristic radii 
of line formation. For Mgxi, Neix, Ovn, and Nvi, this 
resulted in values of 4, 5, 12.5, and 34 stellar radii, respec- 
tively. 

The hotter and higher ionized ions are on average 
formed deeper in the wind, closer to the stellar sur- 
face. This is in agreement with values established by 
Waldron & Cassinelli (2001), based on HETGS obser- 
vations. Applying the standard beta velocity law v(r) = 

Vac (l — T^)^ j with (3 = 0.7(.l) (Grocncwcgen & Lamers 
1989) we obtain at a distance of 4i?„ (where the hot- 
ter plasma of Mgxi is located) v(r) — 0.82 x = 
1540 kms -1 . This velocity is far higher than the veloc- 
ity jump necessary to produce the high temperatures. 

Waldron and Cassinelli (2001) found that line diagnos- 
tics for Si xni indicate that this line emission forms very 
close to the stellar surface, where the velocity is too small 
to produce the shock jump required for the observed ion- 
ization level. The analysis of 17 O stars by Waldron & 
Cassinelli (2007) has led to a strengthening of this con- 
clusion. However, due to the lower resolution of the RGS 
relative to that of the HETGS on Chandra and the weakness 
of the relevant lines in our spectrum, we cannot investigate 
the formation radius of Sixni . These limitations there- 
fore preclude us from confirming or denying the presence 
of this effect in the X-ray emission of £ Ori. 

We conclude that our results are consistent with a 
model of shock fragments that are embedded in the wind 
and are expanding outward with the wind. The FeXVll 
line ratios indicate that the individual source regions as- 
sociated with the shock fragments are optically thin. The 
lack of variability of the X-rays over the years implies a 
large number of sources embedded in the wind. Ions are 
found to be formed and emitting both from far out and 
deep in the wind. The symmetry of the line profiles and the 
slight blue-shift are consistent with the "porosity model" 
(Feldmeier et al. 2003, Oskinova et al. 2006, Owocki & 
Cohen 2006 and Fig. 4 therein), though we find a wider 
range in line formation radii than the w 1.5 i?* used in 
their scenario. 
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parameter 


C Ori Aa 


ref. 


C Ori Ab 

ref. 


COri B 


ref. 


COri C 


ref. 


HR 


1948 






1949 








HD 


37742 






37743 








ADS 


4263 A 






4263 B 




4263 C 




d (kpc) 






fl OK + 0.06 a 
U-^'J_o.04 








ESA97 


E B -v (mag) 






0.06 








VB89 


A v (mag) 






0.19 










logiVH (cm 2 ) (ISM) 






20.34 








FP97 


separation (arcsec) 






0.045 HW00 


2.42 


MG98 


57.6 


MG98 


p (yr) 








1509 


MG98 






spectral type 


09.71b 


WF90 


late O V HW00 


B2III 








9d (mas) 


0.48 ± 0.04 


HD74 












V (mag) 


2.03 


VB89 


4 HW00 


4.2 


HW00 


9 


HW00 


B - V (mag) 


-0.21 


VB89 












(B - V) (mag) 
w rad (kms -1 ) 


-0.27 


VB89 












+25 ±5 


BG78 












v sini (kms -1 ) 


123 ±3 


Pe96 












«oo (kms -1 ) 


1885 


KH96 












T* (kK) 


31.5 ±1 


FP97 












M (M ) 


34 


VB89 












logff (cgs) 


3.2 ±0.1 


VB89 












BC (mag) 


-3.16 


VB89 












log L/L 


5.2 b 


FP97 












M (lO'^oyr^ 1 ) 


1.4 6 


FP97 












Y= riHc/riH 


0.1 


VB89 












An 


normal 


BL82 












T Einstein-ipc ( ke y) 


0.46 ±0.08 


CH89 












logLf nstcin - ipc (ergs" 1 ) 


32.1 ±0.24 6 


CH89 












T x ROSAT (keV) 


0.22 


BS96 












logL* osat (ergs- 1 ) 


32.0 h ' c 


BS96 













Notes: 

a: The distance to £Ori is subject to debate. The three bright Belt stars in Orion, 5 Ori (09.5 II±B0.5 III), eOri (B0 la), and 
COri, form the high-luminosity end of the OriOBlb subgroup, which has an estimated age between 7Myr (Blaauw 1991) and 
1.7 Myr (Brown et al. 1994). The distance to OriOBlb had been estimated to be 360±70pc based on photometry (Brown et 
al. 1994), or 473±33pc based on Hipparcos parallax and proper motion measurements (de Zeeuw et al. 1999). We adopt here 
the direct Hipparcos parallax measurement of C Ori, corresponding to the distance as given above. 
b: Paramater scaled to d = 0.25 kpc. 

c: Berghofer & Schmitt (1994) report that ( Ori showed no evidence for variability over a three year time span, except for a two 
day period in September 23-25, 1992, when an increase of ~ 15% in the X-ray count rate was measured. 
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